Finite Temperature Dynamical Structure Factor of Alternating Heisenberg Chains 
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We develop a low-temperature expansion for the finite temperature dynamical structure factor of 
the spin half Heisenberg chain with alternating nearest neighbour exchange in the limit of strong 
alternation of the exchange constants. We determine both the broadening of the low lying triplet 
lines and the contribution of the thermally activated intraband scattering. 
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I. INTRODUCTION 

Much is now known about the physics of quasi-one di- 
mensional Heisenberg anti-ferromagnetic chain materials 
at zero temperature. This understanding has benefited 
from powerful analytical techniques (see e.g.^i^i^ and 
references therein) as well as highly accurate experiments 
performed at temperatures much smaller than the rele- 
vant exchange constants^. As such the regime in which 
these materials are understood is dominated by quantum 
fluctuations. In contrast, far less is known about their fi- 
nite temperature behaviour where there is an interplay 
of quantum and thermal fluct uat ions^'^'^'^ i ^ o i ^ ^ . 

Though the chemistry of these materials can be ex- 
ceedingly complicated, a wide range of their properties, 
up to energy scales set by their exchange constants, are 
well described by simple lattice models of the form 



H — ^ Jij Si ■ Sj 



(1) 



Here Jy is the exchange constant between spins at sites 
i and j. 

Perhaps the most encompassing probe of the proper- 
ties of spin chains are inelastic neutron scattering exper- 
iments. Such experiments provide detailed information 
on the chains' excitationsiS. In theoretical terms, these 
experiments specifically yield the spin dynamical struc- 
ture factor 



5"^(c., Q) = -i_^-_Im [x"^(c., Q)] , 



TT 1 - e-f^'^ 



x"^(c.,Q) = - 







(2) 
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Here /3 — l/ksT and a, 7 = x,y,z and the brackets 
imply a thermal expectation. We have used the Matsub- 
ara formalism, with imaginary times r and frequencies 
ujn- The dynamical susceptibilities x'^^'i'^^Q) tsike their 
simplest form at zero temperature, where they provide 
direct information on the energies and lifetimes of the 
spin excitations. At finite temperature the x"^{^tQ) 
become more complex functions, determined now by a 



competition between thermal and quantum effects. The 
value of the exchange coupling, J, determines the ex- 
tent of the temperature's role. If J ^ T, there exists 
an appreciable range of energies, u, over which inelastic 
neutron scattering experiments effectively probe the zero 
temperature form of x^'^i^jQ)- However if J '--^ T or 
one is interested in energies, uj, on the order of T, the 
effects of temperature must be taken into account when 
calculating x"^{^^, Q)- 

In this article, we are interested in calculating the dy- 
namical structure factor at low but finite temperature 
for a class of dimerized spin-1/2 chain materials. These 
materials are well described by the Hamiltonian 



N/2-1 

E 

a=0 



( JS2a ■ S2a+1 + «/'S2a+l ' S 



2a+2) 



(3) 



Here J is the exchange coupling through which pairs of 
neighbouring S = 1/2 spins form singlet dimers while 
J' = a J gives the strength of the interdimer interaction. 
We will be interested in the case a <C 1. The lowest 
lying excitations of this model, which we will refer to 
as magnons, are characterized by a gap. A, which at ze- 
roth order in a represents the cost of breaking one dimer. 
With small but finite a, the magnons disperse according 
to 

Ck = J — cos(fcc?), 

where d is the interdimer distance. A particular re- 
alization of this material is Cu(N03)2 • 2.5H2O where 
J = 5.22K and a = 0.27-i^. While this material has 
been studied by inelastic neutron scattering at 300mK— , 
a temperature far smaller than A = AAK , more recent 
experiments have been performedi^ at temperatures on 
the same order as the gap where thermal effects on the 
form of x°'^{^t(1) cannot be ignored. Previously much 
theoretical work has concentrated on the zero temper- 
ature limit, establishing various properties of the spec- 
trum of ([3]), including multi-particle continua and bound 
stateiM^'ia^'21^^ as well as spectral weights^i^ 
and the dynamical structure facto r^^'^^i^^ . However at fi- 
nite temperature, the most pertinent features to explore 
are the broadening of the single particle modes, due to in- 
teraction with the thermally populated background, and 
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the low frequency response arising from intraband tran- 
sitions. 

By virtue of the spin rotational invariance of the 
Hamiltonian ([3|) all off-diagonal elements of the dynami- 
cal susceptibility vanish and the three diagonal elements 
are the same, i.e. x^'^i^iQ) = X^^i^^Q) = X^^i^^Q)- 
In what follows we will therefore only consider x^^- 

Accounting for the effects of temperature in the dy- 
namical susceptibilities makes this problem particularly 
challenging. To see why we consider a Lehmann expan- 
sion of the spin response function in terms of the eigen- 
states of the model, {\l)}. Defining C{t,x) = {S^{t)S^) 



where x — Rj — Rk, this expansion takes the form, 



C(r,x) 



(4) 



The double sum in this representation for C (r, x) arises 
on the one hand from the Boltzmann sum, '^ie~^^^ , 
where Ei is the energy of eigenstate \l), and on the other 
hand from an insertion of a resolution of the identity 
between the two operators, Sj{T) and S^. This expan- 
sion renders the task of finding C(t, x) into a matter of 
computing individual matrix elements {l\Sj{T)\m). At 
zero temperature, this computation is simplified on two 
counts:'^- (i) the sums over eigenstates in Eqn. ^ reduce 
to a single sum; and (ii) the matrix elements needed are 
of a single type, namely those connecting the ground with 
various excited states. At finite temperature however, we 
must deal both with the double sum in Eqn. ^ and the 
matrix elements in their full generality. 

To make this task tractable, we exploit the fact that 
the spin chain material is gapped. On a qualitative level 
the excitations can be divided according to the number, 
n, of magnons they contain. The energy of an excitation 
with n magnons is then at least nA. This notion is im- 
precise because magnon number is not a good quantum 
number, nonetheless at small a it can serve as a rough 
guide to the excitations' energies. In turn, provided the 
temperature does not exceed the gap. A, the contribu- 
tion to the sum in Eqn Q of excitations containing large 
numbers of magnons will be exponentially suppressed by 
the Boltzmann factor, e~^^'. In such a case, we thus 
need only to consider excitations in the Boltzmann sum, 
e~^^' , involving only a few magnons. Concomitantly, 
provided we are interested in determining lmx^^{uJ,Q) 
at energies not far in excess of the gap, we can similarly 
restrict the sum in Eqn.(|3]), X^m' ^-rising from the reso- 
lution of the identity. 

The evaluation of x^^I^tQ) i^i however, more deli- 
cate than the above implies. When evaluating the lead- 
ing terms in the Fourier transform of the spectral repre- 
sentation (ID), one finds divergences when the frequency 
approaches the magnon dispersion. Such divergences 
are expected, since the spectral sum still contains the 
T = result, which is a delta- function at the position 
of the single-magnon dispersion. On physical grounds 
the single-magnon line is expected to broaden at T > 0. 



Analytically this is achieved by carrying out a resumma- 
tion on the divergences of the higher order terms in the 
expansion. Specifically, the sum in Eqn. ^ can be reor- 
ganized according to a Dyson-like equatioujii where we 
write x^^i'-^i Q) in the form. 



DiLu,Q) 



1 - D{uj,Q)T.{uj,Q)' 



(5) 



Here D{uj,q) can be thought of as the propagator for 
non-interacting magnons and T,{uj, q) is the magnon self- 
energy. The key is to match the perturbative expansion 

of e 



X'^iuj, Q) = D{uj, Q) + D^{u, Q)^{uj, Q) + 



(6) 



to the spectral representation of x^^(^7 Q), which is given 
in terms of the Fourier transform of ([4]). In this way we 
obtain a controlled low temperature expansion of the self- 
energy, q), in lieu of x'^^'i^-, Q)- 

This approach has been used successfully in the study 
of finite temperature dynamical correlation functions in 
gapped one dimensional quantum antiferromagnets with 
continuum integrable field theoretic representations^i. 
There the matrix elements, (Z|5|(0)|m), were com- 
puted exactly via analyticity constraints coming from 
integrability.i^^ However, the model of the dimerized 
spin chain, Eqn. ([3]), is not exactly solvable. But be- 
cause a is small, we can compute the necessary matrix 
elements perturbatively in a. 

The T > dynamical susceptibility has been stud- 
ied previously using exact diagonalization of finite length 
chains We believe our approach provides a useful com- 
plement to this work. The numerical approach yields 
results for all a and is not restricted to small tempera- 
tures. However, the system size that can be studied is 
quite small. We, on the other hand, must proceed per- 
turbatively in a and are restricted to low temperatures, 
but our calculations do not suffer from finite-size effects. 
Moreover the nature of low-lying excitations is more ap- 
parent and we can identify the specific processes that 
give rise to the various finite-temperature effects in the 
structure factor. 

With this in mind, a specific feature that we focus 
upon in our analysis is the presence of temperature in- 
duced neutron scattering intensity at low frequencies 
much smaller than the zero temperature gap. The origin 
of this intensity is intraband scattering. The analogous 
phenomenon in Ising-like antiferromagnetic spin chains 
was first pointed out by J. Villain^^ and was first observed 
in the anisotropic spin chain material CsCoBr3<^i2£ For 
Ising antiferromagnets, the relevant excitations are do- 
main walls in the anti-ferromagnetic order. In contrast, 
in the dimer model the relevant excitations correspond 
to low lying magnons. In both cases these excitations 
experience intraband transitions. 

An outline of the paper is as follows. In Section II, 
using first order degenerate perturbation theory, we de- 
termine the dimer model's low lying spectrum and the 
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corresponding matrix elements. In Section III we discuss 
in detail how to use these ingredients to compute the sus- 
ceptibility, x^^I^tQ)- In particular, we explain the use 
of the resummation implied by the Dyson-like equation. 
In the final part of the paper. Section IV, we present the 
actual results for x^^i^i Q)- 



II. GROUND STATE AND EXCITED STATES 
OF WEAKLY COUPLED DIMERS 

Our starting point is the Hamiltonian given in Eqn. ^ 
with an even number of sites, TV, and periodic boundary 
conditions. For small a = J' j J ■, we split the Hamilto- 
nian into a solvable part proportional to J, i?Oi and a 
perturbation iJ' (proportional to J'): 

^72-1 

-ffdimor = {J^la ' S2a+1 + ^'§20+! ' §20+2) 



a=0 



(7) 



We first make some remarks about the J' = case, in 
which the spins decouple into pairs on the bonds J. The 
ground state of is unique and is given by 



Af/2-l 

io>= n 



(8) 



0=0 



|0>. = ^(lT>2aU>2a+l-U>2alT>2a+l) (9) 

SO that |0)q is a singlet between sites and 1. We take 
the associated eigenvalue, Eq — —3NJ/8, as the zero of 
energy. 



A. Excitations 

Excitations are formed by breaking singlets to create 
triplets. The spectrum of Hq then consists of degenerate 
levels at energies nJ relative to the ground state, where 
n < N/2 is the number of triplets. These excitations are 
dispersionless hardcore bosons. When the perturbation 
H' is applied the degeneracies are removed and coherent 
single-particle excitations with dispersion relation are 
formed. These magnons are not free, but interact with 
each other through the perturbation H' in addition to 
being subject to the hard-core constraint. The first ex- 
cited state consists of N/2 — 1 singlets and one triplet, 
leading to a total spin Stot — 1- 

We define da{m) as the operator that breaks a dimer 
between sites 2a and 2a 1, creating a state with z- 
component of spin m. For example the explicit form of 
one of these operators is 



rfa(0)=4(0) = 25|. 



though it is important to realise that d(±l) 7^ d^il). 
A translationally invariant state is formed by taking the 
Fourier transform: 



|p, m) 



7V/2-1 



e'^P-da{m)\0). (11) 



a=0 



Here the factor of two in the exponential accounts for the 
interdimer distance. Periodic boundary conditions lead 
to the quantization condition 



JpN _ ^ 



SO that 



P=—, n = 0,l,2,...,7V/2-l. 

Strictly at the point a = these single particle exci- 
tations are iV/2-fold degenerate with a flat dispersion 
tp — J . A finite value of a causes magnons to 'hop'. 
To first order in a the dispersion is 



J' 

Cp = J - y C0s(2p), 



resulting in a gap 



(12) 



(13) 



(10) 



Using translational and spin-rotational invariance, we 
can express two-magnon states (to lowest order in a) in 
the form 

Ar/2-la-l 

|pi,P2,S',r7l) =7Vs(pi,P2) ^ ^V'ffc(Pl,P2) 

a=l f)=0 

x$frio). (14) 

Here S — 0, 1,2 and the normalisation M will, in gen- 
eral, be dependent on the linear and angular momenta. 
Explicit expressions for the spin part $f^™ are given in 
appendix]^ The wavefunction is given by 

^ab{Pi,P2) = e2'(f^''+^'^'') + <p^e2'(^'^''+P^''). (15) 

Embodied in the non-trivial relative phase, ^p^pj , is the 
magnon-magnon interaction. For a = the form of A 
is unspecified because the magnons cannot hop onto the 
same site and as such do not interact. To lowest order in 
a the correct basis in degenerate perturbation theory is 
given by requiring 

V2H' |pi,p2,S', m) 
J' 

= -y[cos(2pi) + cos(2p2)] \pi,P2,S,m) (16) 

where V2 is the projection operator onto the two-particle 
states. When the triplets in the sum given in (fTil) are 
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well separated (|o — 6| > 1) the condition ((T6l) is triv- 
ially satisfied, independently of A. When the triplets are 



(17a) 
(17b) 
(17c) 



neighbouring 


{\a-b\ 




1) we find 




/I" 

PlP2 


1 + 


e 


-^*tPl+P2j 


2e-2ip2 


1 + 


e 


-2i(pi+P2) _ 


2e-2ipi 


P1P2 


1 + 


e 


-2i(pi+p2) _ 


g-2ip2 


1 + e 


-2i(pi+p2) _ 


e-2ipi ' 


A^ 

PlP2 


1 + 


e 


-2i(pi+p2) _|_ 


g-2ip2 


1 + 


e 


-2i(pi+p2) _|_ 


g-2ipi ■ 



The magnons therefore experience both an infinite onsite 
repulsion and a nearest neighbour momentum and spin 
dependent interaction. Periodic boundary conditions and 
the restriction on the sums lead to the conditions 



<P2 = and <^ 



„ip2JV _ 



(18) 



with the implication that the quantisation of the two- 
particle momenta depends on the total spin, S. These 
are in fact the Bethe Ansatz equations for the spin-^ 
XXZ chain where the sectors S = 0,1,2 correspond to 
anisotropics A = 1,+^,— ^ respectively (and antiperi- 
odic boundary conditions for S = 1). Solving these equa- 
tions for finite N to find pi and p2 is a numerical task 
which we outline in the appendix. We note here how- 
ever that the solutions may be complex, leading to bound 
states, in agreement with results in the literaturei^iiiii^. 
For notational convenience, we define a phase shift by 



'^PlP2 ^ 9 1^ (^PlP2)- 



(19) 



For real pi , p2 the normalization of a two-particle state 
is given by 



N 

2 \ 2 



A/s(pi,P2) = 
N cos{2SlpJ - cos(2pi - 2p2 - 2S^^j,J 



-1/2 



(20) 



2 1 - cos(2pi - 2p2) 

We note that two-magnon states have the symmetry 

\pi,P2,S,m)^e~^<-2 |p2,pi,5,m). (21) 



To avoid an overcomplete basis we make the restriction 
Pi > P2- 



involving at most two magnons. Due to the isotropy 
in spin space of i?dimcr we require only the com- 
ponent of the dynamical susceptibility. In addition the 
Hamiltonian is invariant under translations by two sites. 
When evaluating ([2]) it is then sufficient to consider ma- 
trix elements of the form {l\Sj\m) where j = 0,1 and 
I, m correspond to states with zero, one or two magnons. 
Calculating such elements is simple to zeroeth order in 
perturbation theory and the results are summarized in 
Tables U and [ITl Certain matrix elements are not given 
in the tables because they are identically zero. This can 
be seen by taking account of the fact that the operator 
Si commutes with the total z-component of spin, lead- 
ing to the transition selection rule AS^ = 0. In addition 
some elements are zero by inspection of the states given 
in (in. 



We define the functions 



Us{P,Pl,P2) 

sin(p - pi - 



^J^s{pi,P2)exp{~i{6^^j,^+S^)} 



PlP2 



sm{p - p2 



(5^ 

PlP2 



S^) 



sin(p -pi) 



sin(p - P2) 



(22) 



and 



/iV\2 

Vs's{p'i.p'2.PuP2)^ (yj Ms{pi,P2)Ms'{p'i,p'2) 



X exp 


{*(C2-Ci+(^- 


-S' 




-sin(pi 


~P'l- 


P1P2 ' 


P1P2 


(S 


-5')f) 






sin(pi 








sin(p2 


-P'2 + 


6^ - 

P1P2 


s^: , - 
P1P2 


(S 








sin(p2 


-p'2) 






sin(pi 


-P'2- 


S'^ - 
P1P2 


s^: , - 
P1P2 


{S 





sin( 



sin(p2 - P[ + + ^gp. ~{S- ^Of ) 
sin(p2 ~p'i) 



, (23) 



which are useful when calculating matrix elements that 
involve two particle states. 



B. Matrix Elements 



For small a the gap to excitations is of order J 
and states with n magnons are suppressed by a factor 
exp(— /3nJ) in the thermal trace At low tempera- 
tures, f3J 3> 1, we then make the approximation that 
we only have to consider processes where at most two 
magnons couple to the spin operators. Their contribu- 
tion can be determined by considering matrix elements 



III. SPECTRAL REPRESENTATION AND 
RESUMMATION 

Taking the definition of the susceptibility in the Mat- 
subara formalism ^ it is helpful to expand in terms of 
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operators at even and odd sites: 

•^0 Ll'=0 



(24) 

Using translational symmetry and grouping terms ac- 
cording to magnon number, the susceptibility can be 
written as 

r,s— 

Ja N ^ ^ 

•'^ /,i'=0 



7r,7 



LJn^rj — iuj 

(25) 

Here 7^ is a multi-index enumerating all s-particle states, 
E-y^ and P-y^ are the energy and momentum of the excited 
state 1 7s), Z is the partition function and 

^^7.7e =\{lr\S'o\ls)?+e^'^{lr\S'6\ls){ls\Sl\lr) 

+ I(7.|^il7.>l' + e"'''^(7.|5il7.>(7s|^ol7.)- 

(26) 

We have also suppressed the energy and momentum la- 
bels on Crs for notational simplicity. Carrying out the 
Fourier transform wc have 



4 uj + irj + bj^^ — iv-y 



■My^^^.(27) 



TABLE I: Non-zero matrix elements of the interband type for 
S| acting at sites j = 0, 1. 



(0|5||p,m> 



{pi,P2,0,0 
{pi,P2,2,0 

{pi,P2, 1, 1 
{P1,P2, 1, -1 
{pi,P2,2, 1 

(pi,P2,2, -1 



S||p,0) 
5|b,0) 
S| jp, 1) 
5|b,-l) 
5| ip, 1) 
5,-|p,-l) 



(~l)^\/3l^'^o(P'Pi'P2) 

(-l)'^'\/^t^l(P.Pl.P2) 

(-l)'V^'^l(P'Pl'P2) 
(-l)''''\/^t^2(P.Pl.P2) 



TABLE IL Non-zero matrix elements of the intraband type 
for 5*1 acting at sites j = 0, 1. 



{p',m 


S||p,m) 




Jf{Sm,l — Sm,-l) 


{pi,P2,l,0 


S|bi,p2,0,0) 




3^Vlo(pi,P2,Pl,P2) 


{pi,P2,2,0 


S||pi,P2,l,0) 




3^V2l(pi,P2,Pl,P2) 


{p'l,P2, 1, ±1 


S|bl,P2,l,±l> 


Ti 


J^Vll(p'i,p'2,Pl,P2) 


(p'i,pL2,±1 


S|bl,P2,l,±l> 


±1 


^ -^V2l{p'i,p'2,Pl,P2) 


{p'i,pL2,±1 


S||pi,p2,2,±l> 


±1 


V 773-^22 (pi, P2, Pi, P2) 


(p'i,pL2,±2 


S;|pi,p2,2,±2) 




^ -^V22{p'l,p'2,P\,P2) 



The non-vanishing contribution at T 
from Cio -I- Cqi: 



C'oi 



1 — cos(Q) 



1 



' + irj ~ eg LU + if] + €Q 
(l_e-/3^«)l?(c^,Q). 



is obtained 

(28) 
(29) 



1 



D{uj, Q) is then the bare magnon propagator. The re- 
maining Crs terms, up to C12 + C21, are obtained in a 
simple manner using Eqn. (j27p and the matrix elements 
in Tables U and mi 



Cn = 



1 + cos(Q) 
N 



E 



(jj ■ 



- Ep - ep+Q + ir] ' 



(30) 



Xf/|(Q -f-pi -|-P2,P1,P2) 
1 



^ ' S Pl>P2 

g-/3cQ + pi+P2 _ e~/3(fpi+<!p2) 



^ + eQ+pi+P2 - Cpi - ep2 + 
1 



^ - eQ+Pl+P2 + Epi + ep2 + 



(31) 



We reiterate that the allowed values of pi , p2 in the sum 
above depend on S and that the sum is only over those 
momenta that produce unique two particle states. 

We now make some initial remarks about the structure 
of the Crs terms. Firstly, the two site dimer basis leads 
to a Q dependent prefactor that differs between the in- 
terband terms Cr,r+i,Cr+i^r and intraband terms Crr- 
Secondly, the C12 -I- C21 term will diverge with system 
size as N. This divergence is expected and should can- 
cel with terms arising from an expansion of the partition 
function Z. Lastly, the lu dependence in the denomi- 
nators is such that interband terms diverge as a; ^ eg. 
This divergence is reflected in the intraband terms, which 
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diverge for oj ±J'sin(Q)- Higher order terms have 
stronger divergences. This standard behaviour is a con- 
sequence of the essentially nonperturbative nature of the 
finite T magnon lifetime. Perturbation theory is unable 
to capture, order- by-order, the decay enhancing effect of 
particle-particle scattering processes, hence lifetimes re- 
main infinite to all orders. Treating the interaction accu- 
rately and rendering the lifetimes finite requires a sum- 
mation of perturbation theory terms to infinite order. 
We achieve this infinite summation for certain processes 
as follows. Taking into account interactions between the 
magnons by a Dyson-like equation, we conclude that we 
can write 



D{io,Q) 



1-D{lu,Q)^{uj,Q)' 
Expanding this we obtain 

X^'icJ, Q) = [D{u, Q) + D''{uo, Q)^{uj, Q) - 



(32) 



(33) 



On the other hand, from our low temperature expansion 
of the spectral representation we have 



[a 



01 



C 



10 



Cii + Ci2 + C21 



1 + Zi + Z2 + ... 
(1 - e-^^«) D{uo, Q) + Cii + C12 + C21 + .. 
(1 - e-'3^«) D{iu, Q) + Cii + C12 + C21 

- Zi{l-e-'^"^)D{uj,Q) + ... 

D{uj,Q)+[{Cii + Ci2 + C2i) 

- Zi (1 - e-^'») D{u;, Q) - e-^'^D{u, Q) 



(34) 



The contribution to the partition function from the one 
particle states, Zi = B^^e"^'^", must be included to 
cancel the N dependence of the C12 + C21 contribution. 
Comparing the two expansions, we make the identifica- 
tion 

^{lu, Q) ^D-^{u, Q) [(Cii + C12 + C21) 



(35) 



Finally we calculate the quantity of experimental interest 
as 



1 



lim 

»?^o TT 1 — e^P'^ 



Im 



D{lo,Q) 



1-D{uj,Q)J:{lu,Q) 



(36) 



IV. RESULTS AND DISCUSSION 

We now choose J — 1 and J' — 0.1 so that a = 0.1 is 
small, as required by our expansion. We consider temper- 
atures less than the gap, so that the low magnon number 
approximation holds. We then calculate the dynamical 
structure factor at given lu and Q numerically. Sums 
over momenta such as those in Eqs. ([30|l and ([3T|) are 
performed for systems of N/2 dimers. Analytically, the 
standard procedure for evaluating ([2]) is to take the ther- 
modynamic limit, perform the resulting momentum in- 
tegrals, extract the imaginary part of the resulting sus- 
ceptibility and finally take the limit 77 ^ 0. The zero 
temperature result for S^^{uj,Q) then takes the form of 
a delta function at lo — eq. For the purposes of numer- 
ics, it is necessary to stipulate 77 before performing the 
sums. This results in a broadened zero temperature re- 
sult, a Lorentzian peak of width 77. To obtain accurate 
results at finite T, this width must be small in comparison 
with the thermally activated broadening which scales as 
j'^-PJ ^ In contrast, to avoid finite size effects, the sums 
must be evaluated on a grid in wave vector space that is 
fine enough to resolve the Lorentzian. The condition for 
producing accurate numerics is then 



_Qj 2?7 A-K 



(37) 



In principle by increasing the system size, iV, very small 
values of 77 could be used. However as explained above 
this is only necessary at very low temperatures, where 
the broadening is not sufficiently asymmetric to be inter- 
esting. For the range of temperatures we investigate we 
find a suitable value of 77 to be 0.002. When calculating 
the intraband response we typically use systems of size 
N/2 = 600. The interband response is less sensitive to 
finite size effects and sums using a smaller system size, 
A^/2 = 400, are permissible. We find that the effects of 
using a larger number of dimers are negligible. If the 
delta function in Eqn. (P7)) is ever to be satisfied we must 
restrict the external momentum to Q = Airn/N with in- 
teger n. We leave discussion of further issues affecting 
numerical accuracy, in particular bound state solutions 
of Eqns. (jlSp to the appendix. 

We first consider the behaviour of the one-magnon or 
interband response and examine the lineshape in energy 
by fixing the external wavevector Q. Fig. [T] shows the 
dynamical structure factor at Q = tt for a range of tem- 
peratures. At temperatures less than the gap the main 
feature is a peak at w = eg. At temperatures below 
T 0.3J the peak is approximately Lorentzian. The 
maximum response falls rapidly with increasing tempera- 
ture and the peak broadens asymmetrically. In particular 
the peak becomes skewed, with a tail extending towards 
energies to ^ J. 

The degree of asymmetry is large compared to that 
found for a variety of other spin chains using a semiclas- 
sical approach^i^. Instead it is similar to the asymme- 
try found by exact diagonalizationlS. The broadening 
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also resembles that found for two integrable spin chains 
in Ref. [ll[ which uses the same resummation scheme 
as this paper. Most importantly the asymmetry has 
been found to be consistent with recent data on cop- 
per nitrate^. We also point out that the asymmetry 
is primarily a consequence of the two-magnon interband 
terms C12 + C21. Intraband terms, Cm have vanishing 
spectral weight in this region and their lack of influence 
is confirmed by observing that the one-magnon mode is 
unaffected if we neglect them in the resummation. 

The small peak at lower energies seen in Figs. [U [2] 
and [3] is due to processes that involve bound states. On 
the level of the calculation presented in this paper, the 
bound state response is essentially sharp and grows with 
temperature. This behaviour is an artefact of the order 
at which we truncate the perturbative expansion. Taking 
into account higher order corrections to the magnon dis- 
persion would broaden the peak. Similarly, taking into 
account processes involving e.g. a bound state and a 
magnon in the thermal background and a bound state 
and two magnons in the intermediate state would lead to 
a significant weakening of this feature. 

The behaviour in wavevector space is dominated by the 
static structure factor, 1 — cos{Q), arising from the two 
site dimer basis. As a result, the response is maximal for 
Q = (2n+l)7r and disappears altogether at Q = 2mT as in 
Figs. [2]and[3](for integer n). Exactly at Q = (2n-|- l)7r/4 
the susceptibility, , is a symmetric function of u about 
the point lo = J; however the dynamical structure factor, 
5^^, is not symmetric because spectral weight is shifted 
to higher energies by the factor (1 — e~'^'^)~^. 

The direction of asymmetry in the dynamical response 
at w ~ eg can be understood qualitatively in terms of 
a joint density of states for transitions between occupied 
one magnon and unoccupied two magnon states. For 
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FIG. 2: Wavevector dependence of the one-magnon or inter- 
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Q = the response vanishes. Wavevectors between Q — 3tt/5 
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UJ > Q this takes the form 

(38) 



FIG. 1: The one-magnon response at Q = tt for — 400. 
At T = 0.2 J the lineshape is nearly Lorentzian (inset), but 
is increasingly asymmetric as temperature approaches the en- 
ergy gap. 



Here n(j>) is the thermal occupation number for a one- 
magnon state with momentum p and n{pi,P2) is the 
probability that the two-magnon state characterized by 
momenta pi and p2 is unoccupied. At low temperatures 
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and weak inter-dimer interactions we have approximately 

n{p) e"^"'' , n(]5i,]52) « (1 -e~'^"")(l -e^'^'p^). 

(39) 

For -7r/4 < Q < 7r/4 and 7r/4 < Q < 37r/4 this function 
is skewed towards higher and lower energies respectively. 
On the other hand the specific form of the lineshape is 
dictated by the matrix element M12 and hence by the 
magnon-magnon interaction. The fact that the lineshape 
of the dynamical response at w « eg is skewed towards 
low frequencies in some regions of the Brillouin zone and 
towards high frequencies in others is a consequence of 
the smallness of the ratio of bandwidth to magnon gap . 
This should be contrasted to the findings of Ref. 
for the lineshape in the 0(3) nonlinear sigma model, for 
which the bandwidth is infinite and concomitantly the 
asymmetry was found to always extend towards higher 
energies. 
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FIG. 4: Development of intraband scattering with tempera- 
ture at Q = 7r/2 for N/2 = 600 dimers. 



calculated the next leading contribution C22 using the 
matrix elements in Table and found it to exhibit a 
stronger divergence a,t uj — zLJ' sm{Q). This shows that, 
just as for the interband contributions, a resummation 
needs to be carried out. This is achieved by including 
the intraband scattering contributions in the low tem- 
perature expansion of the self-energy in Eqn. (|32|) . A 
complication that arises in doing so is that because of 
the different prefactors l±cos((5) for intra and interband 
scattering the result of the resummation is reliable only 
at very low temperatures for certain wave vectors. For 
such values of Q higher order terms such as C22 should 
be taken into account. However, a consistent treatment 
of such terms would require the incorporation of inter- 
band contributions involving three- magnon states, which 
is beyond the scope of this work. By including Cn in the 
resummation we remove the square root singularities and 
associated threshold. Instead the response has two finite 
peaks and falls rapidly to zero for \uj\ > J'|sin(Q)|. This 
is physically sensible and is analogous to what is found 
for the intraband scattering in the spin-1/2 Heisenberg- 
Ising chain^^-. In Fig. 2] we show the calculated intraband 
scattering at Q — tt/2. The Q-dependent range in to of 
this scattering compares well with that calculated previ- 
ously by exact diagonalisation of chains of = 16 sites.— 
In that case however, the small system size limited the 
number of available transitions and so the lineshape was 
not representative of the thermodynamic limit. 

In conclusion we have calculated the approximate dy- 
namical structure factor at finite temperature of the al- 
ternating Heisenberg chain in the limits of strong alter- 
nation and low temperature. The method we use has 
previously been applied to integrable spin chains but in 
this case the system is non-integrable. We find that the 
lineshape of the lowest lying one-magnon mode is increas- 
ingly asymmetric with temperature, a direct consequence 
of magnon-magnon interactions. We also establish a pre- 
diction for the low temperature lineshape of the intra- 
band scattering. 



Next we turn to the intraband response. When 
magnons are thermally excited, incident neutrons can 
scatter them within the same band with energy transfers 
small compared to the gap. Accordingly at finite tem- 
peratures there is a spin response at energies w ~ 0. To 
lowest order the contribution to the intraband response 
is given by 



_ 1 + cos(Q) e-^('^-"/^) 

~ 2^ y/{J' sin(g))2~^ 

X cosh cot(Q)v/(J' sin(Q))2 -t^2^ , (40) 

valid for \uj\ < J'\sm{Q)\. This response is bounded by 
inverse square root singularities and has an overall mag- 
nitude that grows with temperature as e~^'-' . We have 
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APPENDIX A: TWO-MAGNON STATES 

The spin part of the two-magnon states is of the form 

{mi,m2} 



9 



The c's are Clebsch-Gordan coefficients but the exphcit 
expressions are given below for convenience: 
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1 r 

1 r 
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diiO)di,{0) 

di{l)di,{-l)-di(^l)di,{l) 
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1 r 

71 

1 r 

71 

1 r 

71 
1 r 
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di{l)di,{Q)-di{0)dv{l) 

dmdi,{-i)-di{-i)di,{o) 
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^l? ^di{i)di,{i), 
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= di{-l)di,{-l). 
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(A2b) 

(A2c) 

(A2d) 
(A2e) 
(A2f) 

(A2g) 

(A2h) 
(A2i) 



APPENDIX B: QUANTIZATION OF THE 
TWO-MAGNON MOMENTA 



In order to carry out momentum sums over two- 
particle states on a finite lattice we require knowledge 
of the allowed values of the momenta, pi and p2 in each 
sector, S. In practice, this means we must solve equa- 
tions (Ull) and (dH) numerically to find the (A^/2 - l)iV/4 
pairs {pi,P2} allowed by the condition pi > P2- This is a 
problem usually encountered in models solvable by Bethe 
ansata^. 



1. Real Solutions 

We first consider scattering states of two magnons, for 
which pi and p2 are both real. In each spin sector S the 
(XXZ Bethe ansatz) equations ([T5| can be written in the 
form 

In order to enumerate all roots of the coupled equations 
(|Bip we take the logarithm. We choose a branch cut such 
that 



= In + ^(2/2 + 1), (B2) 



for S* = 0, 2 and 

Pi 

P2 



AT V P1P2J 'at ' 



N 

i 

'n 



ln( 



^ 

P2pl/ 



N 



:2h 



(B3) 



for S~l. Here the integers /i 2 have range < /i_2 < 
N/2. We note that Ii > I2 implies that pi > p2 and 
using the indistinguishability of particles we can restrict 
ourselves without loss of generality to the case Ji > l2- 
Using the parametrization (|B2IB3|) it is now a relatively 
straightforward matter to determine real roots corre- 
sponding to pairs of integers Ii > I2 by standard nu- 
merical root finding algorithms. 

There are a number of roots which require special 
treatment. Specifically in the singlet sector, there is a 
class of real roots {pi = 2(/i -t- l)7r/A^, p2 — 0} and in 
the quintet sector there is a solution {pi = 7r/2,p2 = 0} 
(for a system with N/2 even). For these cases the deriva- 
tion of equations (|20l22l23p needs to be revisited. For 
the special solutions in the singlet sector the phase shift 
is zero and one finds (setting pi = q) 



AA'(9,0) 



N\^ 



N 



f N N \ 

U{p, q, 0) = U{q, 0) (2 - -Sp,o - Y<5p-,,o ] 



(B4) 
(B5) 



In the quintet sector the same expression is found, but the 
arguments of the two Kronecker deltas above are never 
satisfied. 

The matrix elements for C22 involving the special so- 
lutions are similarly affected and need to be replaced by 

Vs's{q',0,pi,p2) = -Af{q',0)Afs{pi,P2){^ye<''2 
sin(pi ~q' - S^^pJ 



2cos((5^^pJ 
sin(p2 



sin(pi 

P1P2' 



s\n{p2 - q') 



rN\ 

Vs's{q, 0, g, 0) = U{q', 0)AA(q, 0) [-) 



N\-^[N 

4 

2 



yss{Pi,P2,Pi,P2) = 27V|(pi,P2)(y 
sin(pi 



(B6) 



(B7) 



[N 

1 - 

2 



P2 



Kp2) 



sin(pi - P2) 



.(B8) 



The numerical root finder does not converge for 0{N) 
pairs of integers {/f , /|}. Most of these correspond to 
complex solutions of (jB2IB3p . which are discussed in the 
next subsection. So far we have restricted our discussion 
to real roots with distinct integers Ii 7^/2- The reason 
for this restriction is that Ii = I2 correponds generically 
to pi = P2, which does not yield a valid solution of the 
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Schrodinger equation. However, in analogy to what was 
shown in Ref. [s^], there are additional "good" real so- 
lutions with repeating integers Ii — I2, which have to be 
treated with care. 



2. Complex Solutions 

In addition to real roots there exist complex solutions 
of (|B2IB3|) . These give rise to wave functions that ex- 
hibit an exponential decay with respect to the distance 
between the two magnons and hence correspond to bound 
states. As the equations ^TE\\ are closed under complex 
conjugation, complex roots must come in pairs 



pi ^ X + iy , p2 = X - iy , x,y real. 
Adding the two momenta using (jB2|B3P gives 
2 



(B9) 



Pi 
Pi 



-P2 
P2 



2x 
2x 



TT—m 



1 



^j^iIi+12), 



S = 0,2, 
5"= 1. 



(BIO) 



The real component x is now uniquely defined by the 
integers If, Defining (f>s — S{S + l)/2 — 2 and substi- 
tuting for X, the value of y can be found by solving 



S 

Plpt 



using a Newton-Raphson method. The resulting values of 
y are such that ^p^pj or equivalently e^P^^ is either very 
small or very large. Consequently, in order to demon- 
strate that our solutions satisfy both Bcthc Ansatz equa- 
tions (|18p numerically, the value of pi must be known to 
very high precision. Fortunately the matrix elements are 
less sensitive and for the system sizes we consider it tran- 
spires that to evaluate C12 + C21 accurately 17 significant 
figures of pi are sufficient. The two-particle normaliza- 
tion for complex momenta is given by 



Us{x ^iy,x- ty) = [(-l)^e"~| (^^ 

N cosh(yA^) - cosh(4?/ - Ny)i -1/2 



1 — cosh(42/) 



(B12) 



or equivalently 



and for the matrix elements the modification is 

Us{p.Pi,pI) = Y-^s{x + iy,x-iy)e'i'^ 
{(l + A^^pj)cos(2(p-a;)-|5) 

( cos(2p — 2a;) — cosh(2j/)) ^, 



JxN~yN 



+ (-1)' 



-2ix^2y 



(Bll) 



which is an even function of y, as required. 
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